Lagrangian and Eulerian velocity structure functions in hydrodynamic turbulence 
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The Lagrangian and Eulerian transversal velocity structure functions of fully developed fluid 
turbulence are found basing on the Navier-Stokes equation. The structure functions are shown to 
obey the scaling relations Sn(r) oc r 4 " and S„(l) oc inside the inertial range. The scaling 
exponents £ n and £ n are calculated analytically without using dimensional considerations. This 
result is a significant step toward understanding the nature of hydrodynamic turbulence. The 
obtained values are in a very good agreement with recent numerical and experimental data. 
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Understanding of statistical properties of a fully devel- 
oped turbulence from the Lagrangian and Eulerian points 
of view has been a challenging theoretical and experimen- 
tal problem for many years [1-3]. Even a bridge relation 
between Lagrangian and Eulerian structure function ex- 
ponents is a theoretical problem. Also, there has been no 
theory based on the solution of the Navier-Stokes equa- 
tion up to now. Recent progress in numerical calcula- 
tions [4]- [6] and experiments [7] encourages to develop a 
statistical theory. First successful step in this direction 
was undertaken in [8] -[10]. The aim of this paper is 
to elaborate the theory of Lagrangian and Eulerian tur- 
bulence based on the Navier-Stokes equation. We con- 
sider fully developed turbulence in incompressible fluid 
in non-viscous limit and advance substantially the vor- 
tex filament model proposed previously in [8]. On the 
basis of this model we find scaling exponents of velocity 
structure functions and derive a bridge relation between 
the Lagrangian and Eulerian scaling exponents. 

The statistics of turbulence deals with different kinds 
of structure functions. In what follows we are interested 
in Lagrangian velocity structure functions 



#(r)=<|v(t + T)-v(i)n 
and Eulerian transversal structure functions 
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s„ (0 = 



(v(r + 1) - v(r)) x 



(1) 



(2) 



In the first ones, the velocity difference is taken along the 
trajectory of one particle and the average is taken over 
the ensemble of particles. In the second, the velocity dif- 
ference is taken in different points at the same time, and 
averaged over all pairs of points. Since the turbulence is 
assumed to be stationary, both results must not depend 
on time. We restrict our consideration by the values r, I 
that belong to the inertial range of scales. The usual as- 
sumption is that inside the inertial range the structure 
functions do not depend on geometry of a flow, and have 
a finite limit as the Reynolds number goes to infinity. 



In the previous papers we proposed a theory of vortices 
formation built on the Navier-Stokes equation [9, 10]. It 
allowed us to calculate the Lagrangian velocity structure 
functions of different orders [8, 10]. In this paper we de- 
velop and generalize the theory, analyzing a boundary of 
a vortex filament and introducing a notion of a filament's 
age t*. This allows to calculate the transversal Eulerian 
structure functions and to find a relation between these 
two types of structure functions. 

The basic idea of our theory is that vortex filaments, 
i.e. the elongated regions of high vorticity, make the main 
contribution to structure functions. 1 Thus, not decay of 
eddies but their stretching is the main process responsible 
for the formation of structure functions. We will see that 
studying this stretching allows to obtain 'intermittent' 
scaling exponents practically coinciding with experimen- 
tal and numerical results. 

Below, we recall briefly the basic points of our the- 
ory, then we present without derivation the results of the 
analysis of the basic equation and generalize our previ- 
ous results for any (not only 5-correlated) random pro- 
cess. Then we apply the results to find the structure 
functions. 

Following [9], we decompose the spatial derivative of 
velocity into the sum of symmetric and asymmetric parts: 



Here io k is the vorticity, and 6y is a traceless (because 
of incompressibility) symmetric tensor. Differentiating 
the Euler equation, we then get the equations for all 
the components of u>i and bij. Changing to the quasi- 
Lagrangian frame moving along the trajectory of a fluid 
particle, we find the ordinary differential equations de- 
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[1] Besides, this corresponds to the Kolmogorov's law. Actually, 
from Sv oc i 1 / 3 it follows — ~ ( ?' ^ . In the inertial range 



SI <g Lq, hence ui ^> loq. 
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scribing the evolution of and by along the trajectory: 
bij + ji^iUJj - Lo 2 5ij) + b ik b kj + pij = , 

= b n kU>k 

Here py = V^V^p, p is the pressure. The time derivative 
of the second equation gives 

w„ = -p n k{t)u>k (4) 

along the particle trajectory. This equation is a di- 
rect consequence of the Navier-Stokes equation under the 
limit v — » without any additional assumptions. 

Consider a vortex filament where vorticity is high and 
the characteristic radius R of the filament is much smaller 
than the scale of the largest vortices L comparable to the 
size of the flow: 

R ~ (Vw/w) -1 <C L (5) 

In general, p n k is some complicated integral function of 
u. We have shown that inside a vortex filament p n k is 
independent on the local value of ui, and the equation (4) 
becomes linear. 

We now introduce randomness into the equation. Note 
that in (3) we have 3+6=9 differential equations for 
3+6+6=15 values w,, bij and pij. The condition of in- 
compressibility gives one additional algebraic equation 
bu = 0. Hence, both in (3) and (4) we have five free 
functions of time. So, we treat p n k as random functions 
satisfying 

(Pij(t)Pnk{t')) = D ljnk (t - t') 

Since the flow is assumed to be locally isotropic, 

Dij n k{St) = {6 in 5jk + 5 lk 8 : j n )D{8t) . 

Normalizing time, we choose / D(t)cIt = 1. Let D(St) 
be small at time intervals St larger or comparable to 
the characteristic time T g . Neither Gaussian PDF nor 
(5-correlation is assumed. The only important assump- 
tion is that the processes p^ are stationary. Then in the 
limit f > 1 the statistical moments of vorticity behave 
as 

(M n > oce A "*. 

The limit of small values n (A„ « t^ 1 ) corresponds to a 
5-correlatcd process. In this limit A„ obeys a linear dif- 
ferential equation of no more than n-th order with con- 
stant coefficients. We have calculated the values A„ up 
ton~ 100. They appeared to satisfy approximately the 
relation 

A„ ~ an 4 / 3 (6) 

(For exact values of first score of A„ see [8].) The nonlin- 
ear dependence of A„ (6) is the manifestation of intermit- 
tency: {uj 2n ) > (w") 2 . In the limit of large n (A„ > t^ 1 ) 
we find 

A„ ~ An . 



This means that at large values of n intermittency breaks 
up. This allows to find the structure functions of differ- 
ent orders independently. The detailed discussion of the 
properties of the equation (4) see in [11]. 

To illustrate the main ideas, consider a simple model of 
axially symmetric, non-curved vortex filament with uni- 
form vorticity distribution inside its radius R. It contains 
all the main features of a vortex filament. The difference 
from the general case is neglecting the curvature and pre- 
cession of a vortex filament and the gradient of vorticity 
inside it. We shall see below that the first approxima- 
tion corresponds to the late stages of a vortex filament's 
evolution. 

We seek a solution to the Eulcr equation in the linear 
form 

U0 = uj(t)r , v r = a(t)r , v z = b(t)z (7) 

The corresponding pressure is [8] 

p(r,z,t) = P 1 {t)r 2 /2 + P 2 {t)z 2 /2 

From the Euler equation we then find differential equa- 
tions for w,a and b. Combining them, we obtain [8, 10] 

w = -P 2 (t)w 

This equation corresponds to (4). In this model the 
growth of the moments is easy to understand: the ampli- 
tude of a solution remains roughly constant for positive 
P2 and grows exponentially if Pi is negative. Hence, on 
average the vorticity grows exponentially. This moans 
stretching of the vortex filament. From the Euler equa- 
tion for (7) it follows that v r /r = a = —d)/(2u)), hence 
the radius of a particle's orbit decreases as r(t) oc uj^ 1 ^ 2 . 

We note that the cross radius R (5) of a vortex filament 
restricts the region where vorticity is roughly constant. 
The velocity loR on the boundary of the region must be 
less than the large-scale characteristic velocity pulsations 
y : on energetic reasons, the vorticity stops growing in 
the regions where it becomes comparable. Hence, the 
equation (4) is linear only inside the radius R~ Vq/lo; to 
find the solution outside R, feedback effect on pik should 
be taken into account. As the vorticity increases on av- 
erage as a function of time, the radius of a filament de- 
creases, R ~ 1/u). Since the radii r of particles' orbits 
change as 1/y/uj, we have R/r oc 1/y/uJ. This means that 
the approximations of flat vorticity profile and of non- 
curved axially-symmetric filament improve constantly as 
the filament compresses. The vorticity continues to grow 
exponentially (on average) until the cross radius of the 
filament becomes comparable to the viscous scale. Since 
we consider a stationary flow, vortex filaments with large 
radii appear constantly, and at any time there is a set of 
vortex filaments in different stages of their evolution. 

Now we use this model to analyze the transversal struc- 
ture functions (2). For simplicity, we assume circular 
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orbits of particles in a filament. Then 

v(r) = r x lj , = 1 x w 

for any pair of points in the vortex filament separated by 
1. Averaging over all pairs of points in the filament, we 
have (Svl) = l n uj n . 

Now we must take the average over all vortex filaments. 
For given /, the correlation breaks when I becomes larger 
than the radius R of the filament: 



l> R~ Vq/uj. 



(8) 



Hence, calculating S^(l) we must account only the fila- 
ments with lui < Vq. This restriction for vorticity pro- 
duces the restriction for the filament's age t*, which re- 
sults in the set of conditions: 



> 



l m e 



< vz 



for all to. Since A m /m grows as a function of to and 
A TO /m — > A as to — > oo, we get 



ml.' I In 

A m I 



1, V 
A T 



(9) 



It is natural that t* depends on high-order moments, 
since they are more sensitive to rare events in intermit- 
tent media. Now the structure functions take the form: 



S£(l) ~ l n {uj n )(t*) oc l n e A " u oc 

. A„ 
Cn=n-- 



One can see that the cutting parameter Vb affects the am- 
plitudes of structure functions but not the scaling law. It 
means that the boundary of nonlinearity affects on am- 
plitude of structure functions but not their scaling expo- 
nents. The vortex filaments with ages less than i* have 
size larger then I and contribute also to the pre-exponents 
of structure functions. 

This particular example reproduces the main aspects 
of the general case. In the regions where vorticity is high 
the equation (4) is linear, since pik does not depend on 
lu. The boundary of these "linear" regions is restricted 
by energetic reasons: in the pairs of points where velocity 
difference becomes comparable to Vq the equation (4) be- 
comes nonlinear, and the exponential growth of vorticity 
stops. So, the relative velocity inside the vortex filament 
is restricted by the condition 



Av < V . 



(11) 



The orthogonal component of velocity difference in two 
points separated by I inside a filament is Av = lxw, 

The condition (11) then restricts the age of the filament 
by (9), and averaging over all pairs of points gives (10). 

One can use analogous consideration to find the La- 
grangian structure function exponents. In two near 



points of a particle's trajectory separated by time inter- 
val r the velocities differ roughly by Av ~ r w 2 T, where 
r is the parameter corresponding to momentarily curva- 
ture radius of the trajectory. The linear equation (4) and 
exponential growth of vorticity are only valid for small 
relative velocities. The condition Av < Vo then gives 
t < i* , where 

T m (oj 2m ) = T m e A2mt - < V™ 
for all to. Making use of lim (A m /m) = A < oo, we find 



= V /t 



so the Lagrangian structure function is 

^(r) = (A«")cxr n (c 2 "(^))cxr«" , 
c A 2 „ 



(12) 



In [8] we obtained a similar expression for S%(t). It 
is identical to (12) if we choose A = A4/2 ~ 3.1. The 
derivation in [8] used the assumption that the structure 
function S2 was determined by stationary part of the 
probability density function. In this paper we take the 
feedback effect into account. This allows to avoid the 
assumptions of stationarity. 

We do not consider longitudinal Eulcrian structure 
functions in this paper. To calculate them, one has to 
make a complicated analysis of the stochastic equation 
(10) (3) and to investigate the statistical behavior of bij, in 
addition to u. From (3) one can see that some linear com- 
binations of (that do not contribute to the right-hand 
side of the second equation in (3)) are proportional to 
uj. Then the longitudinal scaling exponents may coincide 
with the transverse ones. 

Basing on relations (10) and (12) one can find the 
bridge relation between Eulerian transverse and La- 
grangian scaling exponents: 



A 2n (n - Cn) = 2A„ (n - £„) 



(13) 



Let us now compare the predictions of our theory with 
the results of recent numerical simulations [6] performed 
for Re\ ~ 600. This is the only DNS that presents both 
Lagrangian and Eulerian transverse scaling exponents. 
The scaling exponents are normalized by £ 2 and (3, re- 
spectively (so-called ESS procedure). Our theory has the 
only adjusting parameter A. We use it to fit the La- 
grangian exponent ^4/^2- Then we calculate all other 
values C„ and £„ from (10) and (12). The results are 
presented in Table 1. One can see that the theory is in 
excellent agreement with the experimental data. 

Fig.l presents the Eulerian scaling exponents (normal- 
ized by ESS) as a function of the order n. The theoretical 
prediction is presented together with experimental data 
taken from [6] and [12]. We see that, for high n, the longi- 
tudinal scaling exponents correspond to the theory better 
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than the transversal ones. This is in accord with opinion 
of the authors of [6], who maintain that the results for 
high-order transversal exponents are less certain. 

Table 1 . Scaling exponents normalized by £ 2 in Lagrangian 
and by £3 in Eulerian case. The results of numerical 
simulations are cited from [6]. 



n 


£„/£ 2 (Lag 


range) 


C„/Ca (Euler) 




DNS 


Theory 


DNS 


Theory 


2 






0.71 ±0.01 


0.72 


4 


1.66 ±0.02 


1.66 


1.26 ±0.01 


1.28 


6 


2.10 ±0.10 


2.14 


1.68 ±0.03 


1.74 


8 


2.33 ±0.17 


2.45 


1.98 ±0.10 


2.13 


10 


2.45 ±0.35 


2.64 


2.25 ±0.15 


2.46 



We stress that, to determine A, we used the fourth- 
order Lagrangian scaling exponent. So the theoretical 
prediction for Eulerian structure functions, which is pre- 
sented in Fig.l, has no adjusting parameters. We also 
note that the concavity of the curve corresponds to in- 
termittency of scaling exponents. In our theory it is a 
direct consequence of intermittency of (u n ). 




FIG. 1: Eulerian relative scaling exponents Cn/C3- Theory (• 
+ line) and experiment: transversal exponents taken from [6] 
(V) and longitudinal exponents from [6] (A) and [12] (Oh 

The expressions (10) and (12) allow not only link the 
Lagrangian and Eulerian structure functions of the same 
order, but also to link the scaling exponents of different 
orders. Namely, using (10) and (12) we express A as a 
function of (, n /(,2 and of Cn/ C3 independently for each 
n. Substituting the scaling exponents obtained in [6], we 
find the values of A (see Table 2). We see that all the 
values A obtained from both Lagrangian and Eulerian 
scaling exponents for all measured n coincide up to the 
experimental errors. 

To summarize, we report in this letter of the first sub- 
stantial advance in derivation of both Lagrangian and 
Eulerian scaling exponents basing on the Navier-Stokcs 
equation. We started from the eq. (4), which is a di- 
rect consequence of the Navier-Stokes equation in quasi- 
Lagrangian reference frame. We treat it as a stochastic 
equation, assuming that is a stationary random pro- 
cess. This assumption is natural, since Pij is independent 



of local vorticity. Another important assumption is that 
validity of the growing solutions of (4) is restricted by 
energy limitation, which means that velocity differences 
must be smaller than a large-scale eddy's velocity. Ana- 
lyzing solutions of (4), we find Lagrangian and Eulerian 
transversal scaling exponents. 

Table 2. Free parameter A calculated from the relative 
scaling exponents Cn/C3 (Euler) and £ n /£2 (Lagrange) 



n 


A 






Lagrange 


Euler 


4 


3.76 ±0.13 


3.8 ±0.4 


6 


3.67 ±0.23 


3.5 ±0.3 


8 


3.57 ±0.26 


3.4 ±0.3 


10 


3.54 ±0.35 


3.4 ±0.3 



We stress that the theory developed in this paper has 
only one fitting parameter. No additional propositions 
were used in calculating both the scaling exponents. Nei- 
ther dimensional nor phenomenological hypotheses were 
used. All the obtained relations are consequences of 
the equation (4) derived directly from the Navier-Stokcs 
equation. Thus, the coincidence between the theoretical 
predictions, experimental results and numerical calcula- 
tions is very good. 
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